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TASK  I - OCEAN  MICROSTRUCTURE 


I.  INTRODUCTION 
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I iolcl  obs cT\'citio ns  licivo  revealed  that  the  ocean's  thermoclinc 

possesses  an  internal  structure  consisting  of  a series  of  thin,  laminar  "sheets" 
Oj.  i'’ign  gradient  separated  by  "steps"  ot  only  moderate  density  gradient.  The 
observed  "microstructure"  is  represented  schematically  in  Fig.  1.  The  steo 
height  II  and  the  lateral  extent  of  the  observed  structure  varies  widely  over 
the  range  of  mean  density  gradient  indicated  in  Table  I.  The  observation  time 
rcfci  s to  tnc  time  spent  taking  soundings  at  various  lateral  locations  whereas 
tl^.e  point  persistence  time  refers  to  the  time  duration  between  soundings  at  a 
given  kateral  location. 

Various  mechanisms  have  been  proposed  to  explain  the  observed 
phenomenon.  These  mechanisms  consist  of  double  diffusion,"^  interleaving,^ 
internal  gravity  wave  breaking  and  internal  wave  steepening.  In  this  report, 
the  concept  of  interval  wave  steepening  is  examined.  We  will  first  determine 
i.ae  =et  of  admissablc  v/o.ve  solutions  and  then  examine  the  nonlinear  coupling. 
Comparison  with  field  data  establishes  the  feasibility  of  the  wave  steepening 
approach. 


Table  I - Microstructure  Data  Base 


Observer 


r-  1 > 

Grad  ( — ^ — ) 
' P ' 


Point 

Lateral  Observation  Persistencie 
Extent  Time  Time 


rki'.ve  ti  Tait 


Simpson 


Woods  b.  V/ilc'i 


O (10"'Vm) 
O (10-^/m) 
O (10"'Vm) 


O (35  m)  O (10  mi)  3 day; 


O (10-20  m)  1-3  km 


3 days 


> 33  hr£ 


Density 


k-,  -j- 


H = step  thickness 


K 


\ 


Mean 


Figure  1 Ocean  Microstructure 


II.  INTERNAL  GRAVITY  WAVE  MODEL 


The  equations  of  motion  for  a two-dimensional  incompressible 
fluid  arc: 

Incompressible:  Dp 


Continuity: 


— L L ^ V 
fix  fi  y 


0 


(2) 


where  x is  the  horizontal  coordinate,  y is  the  vertical  coordinate  (positive 
upwards),  and  u and  v represent  the  horizontal  and  vertical  velocities  re- 
spectively, The  superscript  ~ denotes  a dimensional  variable. 

X momentum: 


fi  u 
fit 
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y momentum: 
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rsd 
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(4) 


Equations  (1)  through  (4)  represent  a model  in  which  the  density  may  be  a 
function  of  temperature  and  salinity  only 


P - P (s,  T ) 


(5) 


where  the  Lagrangian  derivative  of  p. 


s and  T are  all  zero. 
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Tlu;  dcMisity  profile  in  a static  aea  will  be  specified  apriori 


Ap_ 

py 


(f') 


where-  is  the  reference  density  at  y = 0.  The  value  of  J is  expressed 
in  units  of  inverse  length  and  has  values  hehveen  10~  Vjti  in  tlie  summer 
thcrmi)cline  and  10  /ni  in  the;  deejj  sea.  The  influc;nce  of  compressibility 
must  be  included  for  valuers  of  ^ of  the  order  of  5 x 10"*^/m.  Thus,  tlie 
incompressible  model  will  be  limited  to  the  thermocline. 


V/e  expand  the  dependent  variables  into  a static  and  wave  contribu- 
tion. 


p^.  Ps  (y) 

P\v 

(7) 

Pr ' \ y ^ 

Ps  ( y)  + p.„  (x»y»t) 

v\' 

(8) 

V ^ y' 

(9) 

y»  t) 

(10) 

v/here  v/e  have  restricted  the  analysis  to  xcro  mean  shear.  The  only  velocities 
are  due  to  the  wave  contribution. 


Setting  all  wave  contributions  equal  to  r.ero,  w'e  obtain  the  hydrostatic 
relationship  between  p and  p . 


^ Pc 
& y 


P^(y) 
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(11) 


From  ((>)  an  1 (7),  we  specify  p 


Ps  - Pr 


(1?.) 


Substituting  Equations  (7)  - (10)  into  (1)  - (4)  wo  w'ill  cletermaio  the  equations 
of  motion  for  the  wave  contribution.  The  equations  are  non-dimensionaliijed 
sucli  tliat  the  linear  terms  are  0(1).  Tlie  reference  variables  arc: 


ref 


where 


I 


P KJ  k 

H V 


o 


X 


ref 


ref 


k 


1; 


>. 

H 

2rr 


II 


X 

y 

2rr~ 


and 


'ref  " '’p  (perturbation  velocity) 


u r = k V / k 
ref  V p H 


^ P V p /a 
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p = V a k P / k . 
ref  p o V r H 


The  nondim ’nsional  equations  of  motion  for  the  wave  contribution  ?^re  ex- 
pressed as 
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V,  her{‘  the  throe  Jimensionles s parauicters  arc 


6 r- 


II 


(17) 
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P /h 
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and 


(19) 


Microstructuro  data  infers  that  both  6 and  a are  much  less  than  unity 
'vhereas  £ may  be  made  small  by  choice  of  the  perturbation  velocity  Vp . la 
ti»c  limit  of  all  throe  parameters  being  small,  the  lowest  order  solution  is 
given  by  the  arbitrary  sum  of  four  plane  waves. 

Lowest  Order  Solution; 
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v.liei'u  A are  arbitrary  constants. 

o , q 

We  seek  a vave  solution  which  propa{-atos  in  the  x direction  and 
possesses  structure  iji  ihe  y direction.  Such  a solution  is  obtained  from  equation 
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1 1 i (t-x  I y) 

2 " 2 ^ 


sin  y 


^i  (t  - x) 


(24) 


Equation  (24)  is  a superposition  of  two  plane  transverse  waves.  The 
veloc:it'_,  jic rturbatio7i  froiti  one  family  of  waves  will  steepen  the.  other  family. 
This  is  unlilvc  acoustic  .sound  v avits  wdiich  arc  longitudinal  w'aves  and  are  self 
steepening. 

In  order  to  describe  the  steepening  of  the  waves  in  equation  (24),  we 
must  include  tlie  lowest  order  non-linear  terms.  Since  the  data  shows 
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a = 0 (10"^^ ) 


wt  retain  (lie  O (e/6  ) terms  in  equations  (13)  - (16). 
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To  solve  cqvtations  (25)  -(28),  we  seek  a solution  of  tlic  same  for;n 
linear  solution 
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(31) 

+ c 

(32) 

where  ^ and  r,  are  constant  on  their  respective  charitcteristics  but  the  char 
istic  must  be  relocated  as  prescribed  by  the  fluid  velocity. 


- t - X - y b F ( ( , Tl) 


(33) 


Tl  - t - X I y + G ( / . n ) 


(34) 


1/ e first  sc'dv  to  dc'seribe 
in  (';/6) 


( f,n)  and  G ( <^ , rj  as  a power  series  expans 
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as  the 


acter- 


e e 2 

F ( .< , r.)  M 5 ) .^1  ( /.  ^ ) t ( ^ ) .<2  ...  p^i) 

e ^ 

G ( Ti)  = ( - ) '1  J , P)  + ( ~ ) ( /,,  ri). ..  (3C) 


SubsMiution  of  (?.9)  - (36)  into  the  equations  of  motion  (25)  - (28),  wo  obtain 


1 


1 i-q 

T 


(37) 


T) 


1 


(38) 


2 


(39) 
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(40) 


Pjquatioiis  (37)  and  (:>d)  illustrate  the  nature  of  trar.s\  ersc  wave  steepenin  '. 
The  velocity  perturbation  from  the  q waves  is  in  the  direction  of  propagation 
of  the  ^ waves.  Hence,  the  r\  waves  steepen  the  ^ waves  and  vice  versa. 


rhe  density  profile  associated  with  transverse  wave  steepeninp;  is 
ootainefi  from  equation  (29).  For  convenience,  we  set  t-x  = TT/2. 


Pw  = y - I (f  ) sin  2 y + O ( y)  (41) 

Equation  (41)  illustrates  that  steepening  can  occur  only  if  6 < 0 (e).  Since  ~ 
is  a small  jjaranieter,  it  is  seen  that  steepening  can  occur  only  for  sn  ail  6 
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and  this  is  in  agreement  with  observations.  The  physical  reason  for  the 

constraint  on  6 lies  in  the  steepening  mechanism  for  transverse  waves. 

Since  the  T]  waves  steepen  the  ^ waves  and  the  T,  waves  are  oscillatory, 

the  ^ waves  must  steepen  within  the  orscillation  time.  Thus  the  perturbation 

velocity  must  disperse  a wave  its  own  w-avelength  in  time  t . This  is  c'-- 

ref 

pressed  as 


V t . > k 
p ref  — V 


-1 


or 


6 < e 


(4?.) 


To  solve  for  the  characteristic  locations  for  arbitrary  values  of  e/5,  we  will 
expand  F and  G in  higher  harmonics  of  e^  and  c^^ 


F ( /,,  Tl)  = 


X J m,  n 


It  ’ ' ' ' ^ ^ 

m.  n 


i m i n Tj 


(43) 


p»q 
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q'H 


(44) 


Where  and  are  constants  to  be  related  by  the  dynamic  equations. 

This  work  will  be  carried  out  during  the  next  quarter  but  the  collective  rosult 
of  this  quarter's  effort  may  be  compared  to  data  to  demonstrate  that  intetnial 
gravity  wave  steep-ning  is  a promising  answer  to  the  source  of  ocean  micro- 
structure. 
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ni.  COMPARISON  WITJrl  DATA 

While  wc  have  yc-t  to  demonstrate  that  steepening  actually  occui-s, 
we  have  sev.-.ral  constraints  which  must  be  consistent  with  observations  of 

microstructure,  otherwise  steepening  will  never  explain  microstrucLure 
data. 

The  condition  for  steepening  requires  s = 0 (6  ).  This  corresponds 


to 


V = 

p 


g p 


2 TT  X 


H 


or,  for  tltc  data  in  Table  I, 


“ 5 A. 

Vp  = 0 (10  m/sec)  to  0 (lO""*  m/sec) 


and 

^ref  ~ ® cm/sec)  . 

Perturbation  velocities  of  the  order  of  cm/sec  are  not  unlikely  in  the 
ocean.  While  observations  demonstrate  that  6 is  0 (lO"^'),  the  present 

theory  demands  only  that  6 be  small  compared  to  unity  ( 6 < 0.  1). 


Furtlicr  constraints  are  placed  on  the  solution  by  viscosity.  We  must 
insure  that  the  inertia  of  the  wave  exceeds  viscous  damping.  This  is  ex- 
pressed as 


or 


0 t 
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(2n)^  V 

^ g P 
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will 


If  the  microsti-ucture  i:;  formed  by  wave  stcepenini*  X 

~ ° V 

corrcapoiid  to  the  step  height  II  and  must  be  smaller  than  the  maximum 
stable  step  indicated  in  Fip„  2.  The  quantity  will  correspond  to  the 
lateral  extent  of  tlie  n'icrostructuro  and  the  wave  period 


(46) 


is  the  persistence  time  of  the  microstructure.  This  is  due  to  the  fact  tlirt 
each  family  of  waves  is  steepened  by  velocity  perturbations  from  the  opposite 

family  which,  in  turn,  is  oscillatory.  Hence,  tlie  micro  structure  will  have 

the  period  T . 

P 


The  analytic  lesults  are  conapared  to  the  data  of  refei’ences  1 throviph 
3 in  Fif.s.  3 to  5 respectively.  Figure  3 represents  a set  of  conditions  under 
which  steepening  is  possible  when  ^ = 0 (lO'^/m)  and  the  data  of  Howe  and 
Tait  lies  within  that  window.  The  theory  indicates  that  the  structure 
should  be  stationary  over  a time  period  of  weeks.  The  data  shows  th  tt  tlio 
microstructure  was  present  and  invarient  for  the  observation  period  of 
days.  Hence,  there  is  no  contradiction  with  data.  Similar  results  are 
illustrated  in  Figs.  4 and  5.  In  all  throe  cases  the  data  falls  within  the 
"window"  under  which  wave  steepening  is  a possible  explanation  of  micro- 
structure. 


Figure  5 


IV. 


CONC.LUSIOi\S  AND  FUTURE  PLANS 


1 he  linearized  'gravity  wave  solution  has 
and  lime  scales  consistent  with  field  data.  The 


demonstrated  length 
conditions  under 


whicn  an  undamped  wave  can  steepen  into 
also  consistent  with  observ "'lions.  Thus, 


a stable  micro  structure  is 
internal  gravity  wave  steepen- 


api^ears  to  be  a promising  prospect  for  the  source  of  ocean  micro- 


structure and  the  research  during  the  i.ext  quarter  will  be  devoted  to 
demonstrating  that  steepening  will  result  in  the  observed  microstructure 


TASK  II  - LASER  SENSOR  TECHNOLOGY  ASSESSMENT 


A variety  of  laser  techniques  for  active,  remote  sensing  of  the  atmos- 
phere and  ocean  luive  been  proposed.  Some  of  these  techniques  involve  Raman 
scattering,  fluorescence,  and  electron  spin  resonance.  From  discussions  with 
A and  IDA,  the  following  areas  have  been  selected  for  study:  water 
teniperature  measurement  and  effluent  detection  using  Raman  Spectroscopy 
and  magnetic  anontaly  detection  using  remote  electron  spin  resonance  (ESJl) 
spectroscopy. 

The  general  detection  equations  can  be  derived  for  each  of  the  above 
techniques.  Also,  there  are  many  sources  of  measurement  degradation  which 
must  be  considered  in  any  specific  scenario.  It  is  the  purpose  of  this  study  to 
try  to  establish  and  document  the  physical  limitations  of  the  detection  processes. 

The  first  area  to  bo  studied  is  Raman  spectroscopy.  It  is  well  known 
that  Raman  spectroscopy  is  a very  selective,  yet  insensitive,  technique  for 
remote  detection  of  gaseous  or  liquid  molecules  and  atoms.  The  general  principle 
in  the  Raman  scattering  process  is  that  the  scattered  light  is  wavelength  shifted 
by  an  amount  dependent  on  the  atomic  or  molecular  structure  of  the  scatterer. 
Thus,  identification  of  the  scatterer  is  possible. 

A numerical  example  will  be  discussed  to  illustrate  a typical  example. 

Let  the  remote  system  be  300m  from  the  scattering  center,  and  the  collecting 
mirror  be  35  cm  in  diameter.  Thus,  the  solid  angle  of  the  receiving  system 
is  10  steradians.  A typical  Raman  cross -section  is  10'^^  cm^/ster.  The 
optical  efficiency  of  the  detection  system  will  be  assumed  to  be  10“^  for  ocean 
surface  detection.  The  attenuation  of  the  laser  beam  and  scattered  beam  must  be 
included  for  sub-surface  detection.  Using  the  detection  equation,  8 the  probability 
of  receiving  a Raman  scattered  photon  per  molecule  for  each  photon  transmitted 
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To  obtain  the  total  probability,  the  number  10"  is  multiplied  by  the 
number  of  scattering  molecules  in  the  resolution  volume.  Taking  an  example 
of  aim  resolution  volume  with  10l6  scattering  molecules /cc,  the  total 
probability  is  10  . Thus,  per  joule  of  green  light,  approximately  100  Ram.an 

Scattered  photons  would  be  received.  As  a reference,  an  average  power  of 
1 watt  at  5030  ^ is  a reasonable  bench  mark  number  for  today's  commercial 
laser  technology.  However,  visible  lasers  with  about  100  watts  of  average 
power  are  under  development  at  several  laboratories.  ^ Maps  will  be  provided 
in  the  final  report  analyzing  the  broad  range  of  possible  systems  parameters 
applicable  for  effluent  detection, 

Raman  spectroscopy  has  been  shown  to  be  capable  of  measuring  water 
temperature,  because  the  structure  of  liquid  water  is  teinperature  dependent. 
r..vperimental  work  and  analysis  by  Chang  and  Young  have  demonstrated  the 
feasibility'  of  this  concept  in  the  laboratory.  However,  field  systems  appear 
to  be  impractical  based  on  todays'  visible  laser  te«:hnology.  Ways  to  enhance 
the  Raman  scattered  signal  are  desired  to  reduce  the  required  average  laser  power. 

Stimulated  Raman  scattering  (SRS)  and  resonant  Raman  scattering  are 
b.-vo  proven  techniques  to  obtain  much  higher  scattering  efficiencies  under  some 
conditions.  It  is  part  of  our  program  plan  to  assess  the  applicability  of  these 
techniques. 

Recent  work  by  Dheer  et  al  has  examined  the  stimulated  Raman  spectra 

12  2 

of  H^O.  Radiation  levels  of  » 100  Mw/cm  appear  to  be  necessary  to  exeeed 
the  threshold  for  stimulated  Raman  scattering.  The  need  for  high  flux  levels  puts 
greater  demands  on  a system  than  a spontaneous  Raman  system  because  of  larger 
peak  power  requirements.  Besides  the  high  flux  level  requirements,  the  SRS 
spectra  is  narrower  than  the  spectra  of  spontaneous  Raman  scattering  because  of 
a gain  phenomena  associated  with  SRS.  The  change  in  spectra  may  influence  the 
sensiii'.-ity  of  the  temperature  measurement  process  and  needs  to  be  evaluator!. 

Our  preliminary  conclusion  is  that  a remote  field  system  utilizing  stimulated 
Raman  scattering  to  measure  water  temperature  may  not  provide  any  advantage  over 
tne  use  of  ordinary  Raman  scattering  techniques.  After  further  study  a more 
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detailed  conclusion  based  on  physical  limitations  will  be  made  in  the  final 
report. 

Daring  the  next  quarterly  period  work  will  be  initiated  on  investigating 
the  use  of  SRS  and  Resonance  Raman  Enhancement  to  determine  the  funda- 
mental limitations  of  this  technology  to  the  detection  of  various  species  such 
as  oil,  H^,  and  other  species  in  sea  water.  In  addition,  the  feasibility  of  ESR 
spectroscopy  to  remotely  detect  small  changes  in  the  earth's  magne  ic  field 
will  be  studied.  The  conclusion  of  both  of  these  analyses  will  be  included  in  the 
Contrac . Final  Report. 
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